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Weak solutions to lubrication equations in the presence 

of strong slippage 
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Abstract 

The existence of global weak solutions is proved for one-dimensional lubrication 
models that describe the dewetting process of nanoscopic thin polymer films on hy- 
drophobyzed substrates and take account of large slippage at the polymer-substrate 
interface. The convergence of these solutions as either the Reynolds number or the 
capillarity goes to zero, as well as their limiting behaviour as the slip length goes to 



^ • zero or infinity are investigated. 



1 Introduction. 



During the last thirty years lubrication theory was successfully applied to modeling of 
)0 • dewetting processes in micro and nanoscopic liquid films on a solid polymer substrates see 
e.g. Oron et al. [H, Miinch et al. [11 and references therein. The influence of intermolecular 
interactions that are typically due to the competition between the long-range attractive 
van der Waals and short-range Born repulsive intermolecular forces play an important role 
in such processes, see Oron et al. [H, de Gennes [sj. 

Besides intermolecular forces and surface tension at the free surface of the film it has 
^ ■ been shown by Fetzer et al. (2] that the dewetting of polymer films on hydrophobic sub- 
strates also involves such boundary effect as slippage on a solid substrate. The measure of 
slip is a so-called slip length, which is defined as an extrapolated distance relative to the 
wall where the tangential velocity component of the liquid vanishes. Recently, it has been 
shown experimentally and theoretically that the early stages of the dewetting process and 
the evolving morphology depend markedly on the magnitude of the effective slip length. 
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which can be of the size of the height of the hquid film or even larger for nanoscale systems, 
see e.g. Fetzer et al. jsf, Miinch and Wagner j^. Recently in Miinch et al. [21], Kargupta 
et al. j7| closed-form one-dimensional lubrication equations over a wide range of slip lengths 
were derived from the underlying equations for conservation of mass and momentum, to- 
gether with boundary conditions for the tangential and normal stresses, as well as the 
kinematic condition at the free boundary, impermeability and Navier-slip condition at the 
liquid-solid interface. Asymptotic arguments, based on the magnitude of the slip length 
show that within a lubrication scaling there are two distinguished limits, see Miinch et al. 

i. 

These are the well-known weak-slip model 

dth = -d, ((/i^ + b h^)d, {ad^^h - U{h)) ) (1.1) 
with b denoting the slip length parameter, and the strong-slip model 

Re {dt{hu) + d^^hu")) = Ad^{u{h)d,u) + hd^iad^^h - U{h)) - - (1.2a) 

dth = -d,{hu), (1.2b) 

respectively. Here, u{x,t) and h{x,t) denote the average velocity in the lateral direction 
and the height profile for the free surface, respectively. The positive slip length parameters 
b and /3 are related by orders of magnitude via b ~ ?7^/3, where the (small) parameter rj, 
< 1] 1, refers to the vertical to horizontal scale separation of the thin film. The high 
order of the lubrication equations (11. ip and f ll.2ap - fll.2b|l is a result of the contribution 
from surface tension at the free boundary, reflected by the linearized curvature term adxxh 
with parameter a > 0. A further contribution to the pressure is denoted by Il{h) and 
represents that of the intermolecular forces, namely long-range attractive van der Waals 
and short-range Born repulsive intermolecular forces. A commonly used expression for it 
is given by 

1 a 

Il{h) = ^ - ^ with a > 0. (1.3) 

The terms Re {dt{hu)+dx{hv?)) , with Re > denoting the Reynolds number, and Adx{y{h)dxu) 
in fll.2ap - (11.2bp represent inertial and Trouton viscosity terms, respectively. We assume 
below linear dependence of the viscosity coefficient on the height, namely 

v{h) := uh, u > 0. 

For these choice of the constitutive laws, we first investigate the existence of global weak 
solutions to fll.2al) - (]1.2bp on the interval (0, 1), supplemented with the boundary conditions 

u = at s = 0, 1 , (1.4) 
dxh = at X = 0, 1 , (1.5) 

and initial data (/io,mo) with positive first component ho- Observe that, in particular, the 
boundary conditions (11.40 for u guarantee the conservation of mass 



/ h{x,t) dx = he := / hQ{x) dx , t > 0, (1.6) 
Jo Jo 
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where he is the average of the height profile. 

We next investigate the behaviour of these weak solutions as some of the parameters in the 
model goes to zero or to infinity. More precisely, let us point out that the strong-slip model 
f ll.2ap -f ll.2b|) includes as limiting cases three further lubrication models. One is obtained 



from the strong-slip model in the limit /3 — )■ oo and describes the dynamics of suspended 
free films, see e.g. Brenner and Gueyffier [sj : 

Re (dtihu) + d,{hu^)) = AdMh)d.u) + hd,{ad,,h - U{h)), (1.7a) 
dth = -d^{hu). (1.7b) 

In the second one, we neglect the inertial terms which corresponds to a vanishing Reynolds 
number. For the third limiting case which is derived in Miinch et al. Q the slip length 
parameter /?/ is of order of magnitude lying in between those that lead to the weak-slip 
model fll.ip and the strong-slip model fll.2ap -( ri.2bp . i.e. 6 ^ /3/ <^ /3. The corresponding 
intermediate- slip model is given by 



dth = -d, (h'^d^ {ad^^h - U{h)) ) . (1.^ 



It can be obtained by rescaling time in f ll.ip by b and letting 6 — )■ oo or by rescaling time 
and the horizontal velocity by P in fll.2ap - fll.2bp and taking the limit /3 — )■ 0. We consider 
here the second limit and change variables in fll.2ap - fll.2bp as follows: 



— — u 

X := X , t:=t(3, h := h , and u := —, (1.9) 

/-^ 

where x, i, h, and u denote the old variables. We show that, as /3 — ?■ 0, the rescaled 
solutions to fll.2ap - fll.2bp converge to a solution to 

u = hd,{ad,,h-n{h)), (1.10a) 
dth=-d.,{hu), (1.10b) 

satisfying boundary conditions fll.4l) - fll.5p that is, a solution to the intermediate-slip equa- 
tion ffTSj) . 

Our proofs follow a strategy that has also been employed in other systems arisin g in 
fiuid dynamics and phase transition theories. In fact, Bresch et al. joj and Bresch jloj | 
considered the following Korteweg system 

dt{hu) + dw{hu®u) = div (/i(/i)(Vu + (Vn)'^) + V(A(/i)div n) (1.11a) 

+ ahVAh - VP{h) 
dth = - div (hu) . (1.11b) 

In the case 

a = 0, fj,{h) = fj, = const, \{h) = A = const. 
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the svstem ffTTTaD - fll.llbp gives the compressible Navier-Stokes equations. For the latter, 



considered on a bounded domain with P{h) having a globally Lipschitz continuous deriva- 



tive, existence of a unique strong solution was shown by Solonnikov [Uj in dimensions 
d > 2. Recently, Bresch et al. jol showed existence of a weak solution to fll.llap -f ll.llb|) 
on bounded domain with 

a > 0, fi{h) = uh, \{h) = (1.12) 
and similar requirements on P{h) for dimensions d > 2. In the proof the authors use a 



so called BD entropy equality introduced by them originally in Bresch and Desjardins [12 
for the case a = 0. In the latter case f ll.llap - fll.llb|) coincides with viscous shallow-water 
equations. In a recent review by Bresch [10] on those equations it was shown that the BD 
entropy equality holds for f ll.llap - fll.llbp in dimensions d > 2 if a special relation between 
viscosity coefficients is satisfied, namely 

X{h) = 2{n'{h)h- fi{h)) (1.13) 

that is satisfied for f ll.l2p . 

We also notice that the strong-slip lubrication equation fll.2ap - (]1.2bp has a similar form 



to the one dimensional case of f ll.llap - lll.llbp . Indeed, in one space dimension (11. 1 lap - 
I ll. llbp reduces to 



dt{hu) + d^{hu^) = d^{u{h)d^u) + ohdlh - d^P{h) (1.14a) 
dth = -d^{hu). (1.14b) 



Recently, in Mellet and Vasseur [13[ an analogue of the BD entropy identity was shown 
for fll.14ap - fll.14bp with any smooth viscosity function ulh) and cr = 0. Furthermore, 
they proved existence of strong solutions to fll.14ap -f ll.14bp on the whole real line M with 
(7 = 0, u{h) = uh^ and P{h) = h"', where A; < 1/2 and 7 > 1. 

The main difference in f ll.2ap - fll.2bp compared to f ll.l4p is the special form of the inter- 
molecular pressure (11. 3p which is singular when the height h vanishes and complicates the 
analysis. Nevertheless, it also provides a very useful positive lower bound for h. Equations 
f ll.2ap -f ll.2bp also include an additional slip term which is easily handled for the existence 
theory but plays a crucial role in the limit /3 — )■ (leading to equation (11. 8p ). 

In this paper we start in section |2] with recalling the energy identity associated to 



f ir^ -f lObj) and derive a version of a BD entropy identity. These two estimates provide 



the main ingredients of our existence results. In section [3] we follow the strategy of Bresch 



and Desjardins [IJ] by setting up a higher-order regularized equation for which analogous 
energy and entropy inequalities are satisfied. The a priori estimates allow us to pass to 
the limit in the regularization parameter to obtain global weak solutions to fll.2ap - fll.2bp . 



In section H] we establish existence of global weak solutions for the cases Re = 0, /3 = 00 
and (T = by passing to the respective limits in f ll.2ap - fll.2bp . While the first two cases 



are straightforward, we need to refine the estimates from section [2] in order to enable us to 
deal with the limit a — )■ 0. Finally, in section |5] we establish that solutions of f ll.2ap - fll.2bp 
converge, after the rescaling (11. 9p . in the limit /3 — )■ to a solution of the intermediate-slip 
model flTTg]) . 
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2 A priori estimates. 



In this section we state two relations that are satisfied by classical solutions of fll.2ap - 
(]1.2bp . As it was stated in the introduction both have their counterparts for viscous 
shallow- water equation fll.llap - (]l.llb|) and lubrication equations (II. ip and (11. 8p . For the 



latter they were initially derived by Bernis and Friedman [15|. The energy equality for 



shallow-water equations is known already for several decades, whereas the entropy equality 
was suggested recently by Bresch and Desjardins [l^. For the strong-slip model fll.2ap - 



|1.2bl) the derivation of the energy equality is again standard, see Kitavtsev and Wagner 
IgI . As a new result we derive here the entropy equality for fll.2ap -f ri.2bp following the 



approach of Bresch and Desjardins [l2j. At the end of this section we use the energy and 
entropy equalities to derive a priori estimates in the case a > on classical solutions to 
fll.2ap -f ll.2bp having positive h and satisfying (11.61) . fll.4p -f lL5|) . For consistency we start 
with the energy equahty. 



Lemma 2.1 (Energy equality). For classical solutions of the system U.2a\) - fl.2b\) with 
boundary conditions ( [i.y^[ )-( l773]) the following equality holds 



d 



Re4 + C/(ft) + .M^ 



dx 



h'h\dxu\'^ dx 



^dx, 



(2.1) 



where the potential function U is the indefinite integral of U defined by 



U{h) 



1 a 

+ 



2/i2 3/i3' 



/i > 0. 



(2.2) 



Proof: One way to show (12.11) is to use the fact that 



E{u,h) :-- 



u 



Reh— + U{h)+a 



dx 



(2.3) 



is a Lyapunov functional for the system fll.2al) - f]1.2bp . see Kitavtsev and Wagner [16]. Here 
we give another standard derivation of (12.11) . Multiplying f ll.2ap by u and integrating in x 
we obtain 

Re / \dt{hu)u + dx{hu^)u\ dx — A dx^i^ hdxu)u dx 
Jo Jo 

— I hudx{(jdxxh — 11(h)) dx + / — dx = 0. 
Jo Jo P 

Using now several integrations by parts, fll.2bp and boundary conditions fll.4p -f fT3|) we 
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obtain 



Re 



dt ( /i-y") + Y - hu^d^u 



dx 



1 rl 
2 



+ 4 / iyh{dxu) dx— / dth^ad^xh — Ii{h)) dx + I dx 
'o io io P 



Re 



dt ( h— \ - dx{hu) — - hudx ( y 



dx 



+ 4 / z//i(9^M)2dx+ / dt(a^-^^ + U{h)] dx+ I ^dx 



u 



Re / dt\h—]dx + A ph{dxuYdx+ dt \ a'-^-^ + U{h) ] dx + — dx 



'0 \ ^ / ^0 
From the last expression (12. ip follows. 



/3 



Lemma 2.2. For classical solutions of the system U.2a]) - fl.2b\) with boundary conditions 
( [j.^[ j-( l773]) the following equality holds 

ld_ f^\dj^ 
2dt Jo h 



dx = j dx { — ) dxU h dx. 



h 



(2.4) 



Proof: The statement is verified again using several integrations by parts, equation f ll.2ap 
and boundary conditions fll.4p - fll.5p as follows. 

1 d \dxh\'^ 
2di 



h 



dx 



'dxhdxth. 1 f'dth\dxhW 
dx / — dx 



h 







dx \ ~J^] dx{hu) dx - j udxhdx 
dx I ) dxuhdx. 



dxh 



dx 



dxh\ 

dx 



Let us define a so-called entropy function by ip{h) := 4z/log/i for h > 0. Then we have the 
following lemma. 

Lemma 2.3 (Entropy equality). For classical solutions of the system lll.2a\) - fl.2b\) with 
boundary conditions ( [i./^[ )-( l773]) the following equality holds. 

d 



dt 



h{Reu + dxm? - + Re (^a^ + U{h) 

-Re f ^dx-Aav [ \dxxh\^ dx - Av f Yi' {h)\dxh\^ dx. 
Jo P Jo Jo 



dx 



(2.5) 



Proof: Let us multiply equation fll.2ap by d^ip^h), integrate with respect to x and use an 
integration by parts, f ll.2b|) and fll.5|) : 



AuRe / {dfU + udxu)dxhdx = —IGu'^ / 
Jo Jo 

Jo Ph Jo Jo 



d^u h dx 



{h)\dM^dx. 



(2.6) 



On one hand, a further integration by parts in the left-hand side of (12. 6p . equation (]1.2bp . 
and the energy equality (12. ip give 



4i/Re / {dtu + udxu)dxhdx 
Jo 



4i/Re 
4i/Re 
4i/Re 
Re 



d 
dt 
d 
dt 
d 
dt 



udxhdx— / udxthdx+ / ud^udxhdx 



udxhdx— / dxudx{uh) dx + / udxud^hdx 



udrh dx 



h[dxuY dx 



dt 



Auudxh + Reh^ + U{h) + (T^^ 



dx + 



dx 



On the other hand, using Lemma 12.21 and the definition of ip, one can write the first term 
on the right-hand side of (12. 6p as 

-16u^ j (^~]~^ ^^^^ = j h\dxip{h)\^ dx, 

while it follows from ( ll.2bp and (II. 4p that the second term on the right-hand side of (12. 6 p 
can be transformed into 



-4z/ 



^ udxh 
~W 



dx = —Au 



I3h 



dx + 



1 d 



/3 "^"^ Pdt 



(p{h) dx. 



Substituting finally the last three identities into (12. 6p one obtains 



d_ 

diJo 



Re" 



■hu + -h\dx(p{h)\ + AuReudxh — —(p{h) + Re [ a 
Z p 



\dxh\ 



+ U{h) 



dx 



= -Re f ^dx-Aau [ \dxxh\^ dx - 4u [ U'{h)\dxh\'^ dx. 
Jo P Jo Jo 

Using the definition of (p the last expression can be easily transformed into (12. 5p . 
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Lemma 2.4. For smooth functions {h,u) with a positive first component h, we have 



h\dx(p{h)\^ dx < 



Re 

hiReu + dMh)f dx + Re Eiu, h) H 



where E{u, h) is defined in Ili2.3\) . 
Proof: Using the elementary inequality 



the fact that 



U{h) > -l/{6a^) 

for all h > 0, and the definition (12. 3 p of E one obtains 



h (Re u + dx^p{h))'^ dx > 







1 

> - 

- 2 

1 

> - 

- 2 



h\d^^{h)\'^ dx - 2Re 



h\dx^{h)\'^ dx — Re^ j hu^ dx 

|2~ 



E{u,h) 



dx 



h\d^^{h)\'^ dx - 2Re E{u,h) 



Re 
3^ 



(2.7) 



(2.^ 



from which the statement of the lemma follows. 



We now deduce several bounds from Lemmas 12.11 and |2 . 31 for classical solutions {h,u) to the 
system ( ll.2ap -( ll.2bp with boundary conditions (I1.4p -( 1T3]) and a positive first component. 
To this end, we assume that the initial data and mq for (ll.2ap - flL2bl) satisfy 

ho e H\0, 1), ho{x) > for all x E [0, 1] and Uq E ^^(0, 1). (2.9) 

As already mentioned, a first consequence of (]1.2bp and flL4l) is the conservation of mass 
f lL6p . We next combine (11. 6p . the energy and entropy identities to obtain the following 
bounds when a > 0. 

Proposition 2.5 (a priori estimates). For fixed positive a,Re,P,T and initial data sat- 
isfying ^2. 9\) there exists a positive constant Co > 1 depending only on T , a, v , Re, a , 
E{uq, ho), and Sp{uo, ho) defined in ^2.13\) below such that the following terms are bounded 
by Co in respective norms 



Vh, d^Vh, h-^^\ d^h, VR^Vhu E L°°{0,T; L\0,1)), (2.10) 
d^{h-^/^),d,,h,Vhd,u,^ E L2((0,1) X (0,T)), (2.11) 

and 

Co^ <h{x,t) <Co (2.12) 
for all X E (0, 1) and t E (0, T). The constant Co tends to oo as a ^ 0. 
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Proof: Integrating the energy equality (12.1 p with respect to time over (0, t), t E (0, T), 
and using the inequahty 



ensure that 



Re,/J^uit)\\l + ^\\h-'/'it)\\l + \\y^dMml < Ci:=2E{uo,ho) + 



3a 



2' 



4:u\\^/h{s)d:^u{ 



+ 



u{s) 


2" 


VP 


2 



ds < Ci. 



Therefore, the norms of ^/a^xh, h and ^/RB^/hu in L°°(0, T; -L^(0, 1)) are bounded by a 
constant that depends only on Ci while the norms of y/hd^u and u/^/^ in L'^{{0, 1) x (0, T)) 
are bounded by a constant that depends only on Ci and i/. Using the bound on dxh one 
obtains 



\h{x,t) - h{y,t)\ 



dxh{z, t) dz 



< \x-y\"^ ||W)ll2< 



for all {x,y) G (0, 1) x (0, 1) and t G (0,T). Integrating the above inequality with respect 
to ?/ G (0, 1) and using (11. 6p readily give the upper bound in (l2.12p . To establish the lower 
bound for h in (I2.12p . we combine the L°°(0, T; L^(0, l))-estimates on and \fadji 

just established to obtain a bound on the norm of 1/V^ in L~(0, T; W^^^{^, 1)) smce 



Wdx{h-'i^)\dx = \l 

Jo ^ Jo 



1 \dM 



h3/2 



dx < \\^dxh\\2 \\h-^'^\ 



2- 



Owing to the continuous embedding of 14^^'^ (0, 1) in L°°(0, 1), the positive lower bound in 
f l232|) readily follows. 

Next, let us introduce the functional Spiu^ h) defined by 



Sp{u, h) :-- 



1 



\dM' 



-h^Reu + d^i^^h))^ + -{Avh-^{h)) + Rea'-^-^ + ReU{h) dx. (2.13) 



It follows from the mass conservation (11.61) and the entropy equality (12. 5p that 



Sp{u{t)Mt)) 



lo Jo (y"' + ^^^l^--^l' + ^''n'(/i)|9,/i|'^ rfxrfs = S';3(Mo,M (2.14) 



for any t G (0,T). Since z — log(z) > 1 for all 2; > 0, it follows from the energy equality 
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(1^ . fETH]) . and Lemma El that 

Re 



Sp{u,h) >\! h{Reu + d:,(f{h)f dx + ^ f {Avh - ip{h)) dx - 
2 Jo P Jo 



6a^ 



> ]- f\{Reu + dMh)fdx-p^ (2.15) 

> ^ /" /i|9,¥p(/i)|2t/x-ReE(M,/i)-:^ 
4 Jo 

> ^ /" /i|9,y.(/i)|2rfx-ReE(Mo,/io)-Tp^. 
4 Jo 

Combining this with the previous estimate on y/adxh, f l2.14p . and the fact that 



n'(/i) > 



2a /eV 1 



imphes that 

\ h{t)\dMhm' dx + 1^ (^4au\d.Ms) r + ^''"[^^^y^'' ) dx ds 



U — / \dMs)Vdxds<a^ . . , ^ 

with C2 ■= S^{uo,ho) + ReE{uo,ho) + (Re/3a^). Since h\dxip{h)\'^ = 4:\dxVh\'^, this com- 
pletes the proof of flfrUjl - fimjl . ■ 



Remark 2.6. As in Theorem 1 of Bertozzi et al. [17[ that shows existence of uniform lower 
bound for solutions of fll.Sp many estimates in the proof of Proposition 12.51 heavily rely on 
the positivity of the parameter a in 11 and become unbounded as a — )■ 0. The estimates 
(12.1 op . (12. lip , and (I2.12p are no longer valid in that case. The same remark holds for the 
case (7 = since the estimates depend on the bound for ^/a^xh. On the contrary the proof 
holds without changes for the cases Re = or /3 = cxd except for the fact that then one 
looses estimates on Re\/hu and u/a/^ in (I2.10p - (l2.1ip . respectively. 



3 Existence of weak solutions. 



We first define a weak formulation of the problem (ll.2ap - fll.2bp with boundary conditions 



(ll.4p - (ll.5p . Consider two functions and satisfying (12. 9p . 

Definition 3.1. A pair {h,u) is a global weak solution to U.2a\) - l[l.2b\) with boundary con- 
ditions ([i.^|j-(fTr3]) and initial conditions {hQ,uo) if h and u enjoy the regularity properties 
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stated in li2.10\) - f2.11\) - (2J^) and the following holds 



^0 



hdtil> dxdt + / ho^p{■,0) dx = — 



oo fl 



^0 



hudxip dxdt, 



(3.1) 



Re / / hudt(j) dxdt + Re / /io%0(-, 0) rfx + Re / / hu^dx(j)dxdt 
Jo Jo Jo Jo Jo 



+ 



oo fl 



hdxudx4> dxdt — a 



JO 
oo fl 



oo fl 



dxhdxxh(p dxdt — a 



JO 



Ili{h)dx(f)dxdt — — 

P Jo 



JO 
oo fl 



oo fl 

Jo 



ud) dxdt = 



for all ip G C^{[0, 1] x [0, oo)) and (f) G C^{{0, 1) x [0, oo)), where 

POD 

Ui{h) ■=- Tl\'{T)dT. 
Jh 



hdxxhdx4> dx dt 
(3.2) 

(3.3) 



To show the existence of weak solutions to fll.2ap -( ll.2bp - fll.4p -f ll.5p we construct ap- 
proximating systems similar to those suggested by Bresch and Desjardins 1^ for the Ko- 
rteweg and viscous shallow- water equations fll.llap -f ll.llbp . Although the pressure term 
Il{h) does not need a regularization as in the case of Bresch and Desjardins lj|, one still 
needs to regularize the function h sufficiently in order to control additional higher order 
terms arising in the entropy equality. The approximating systems we take are given by 



Re {dt{heUe) + dxQleUl)) = Adx{vhedxUe) + hedxicrdxxhe - n(/le)) 



(3 



(3.4a) 
(3.4b) 



dtK = - dx (heU,) , 

where £ > is a small parameter. Consider f l3.4ap -f l3.4bp with boundary conditions 

Us = dxxUe = dxK = dike = d% = 0, (x, t) G {0, 1} X (0, T), (3.5) 

and initial conditions 

Us{x,0) = Ui;^o{x) and /^^(x, 0) = /ie_o(2^) > 0, x G (0, 1), (3.6) 

where u^^o and h^fi are smooth functions such that 

Ue,o Uq in L^(0, 1), h^fl Hq in H^{0, 1), and eh^fl in H^{0, 1) as e 0. (3.7) 

Proceeding as in Lemmas 12.11 and 12.31 one derives formulas analogous to (12. ip and (12. 5 p 
for the system f l3.4ap - fl3.4bp with boundary and initial conditions fl3.5p -f l3l6|) . Namely, 
introducing the energy 



E,{us,hs) : = 

Jo 



\dxhe 



u 

Reh,^ + U{h,) + a^ 2 



dx. 
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and the entropy 



Se{Ue,he] 



-he{ReUe + dxip{he)f - -^Vihe 



+ Re 



a- 



+ U{h,) + -\d% 



dx 
dx, 



the corresponding energy and entropy equahties read 



^Es{ue, he) = -A j i'he\dxUe\^ dx - j ^dx-e^ I \dncxUe\^ dx 



(3.8) 



and 



-Re f ^dx-Aau ! \dxxhe\^ dx - 4u [ W {he)\dxhe\'^ dx 
Jo P Jo Jo 

- [ [Re€^\d,,Ue\^ + Au€\d%\^ + 4:ue^dxxUedl\oghe] dx. (3.9) 
Jo 

Given e > equation fl3.4ap is parabohc with respect to u. Arguing as in Bresch and 
Desjardins [ij], the initial-boundary value problem system f l3.4ap - fl3.4bp with fl3.5l) - fl3.6p 
has a unique classical solution at least locally in time. The next proposition establishes a 
uniform lower bound for he and thus guarantees the global in time solvability for fl3.4ap - 
([330). 

Proposition 3.2. For fixed positive e, a, Re, /3 and T > let {he, Ue) he the solution of 
the initial-boundary value problem system \3.Ji.d^ - [3.Ji.b\) with hS. 5\) - [3^) in (0, 1) x (0,T). 
There exists a positive constant C3 > 1 depending only on T , a, v, Re, /3, a , ho, and 
Uq such that, for all sufficiently small e > 0, the following terms are bounded by C3 in 
respective norms 

(3.10) 
(3.11) 



dx {h;'/^) , dxxhs, ^/hedxUe, Ue, ^ed%, sd^xUe G ^^((0, 1) X (0,T)), 



and 



^<he{x,t)<C-i, (x,t) G (0,1) X (0,T). 



(3.12) 



Proof: Throughout the proof C denotes a positive constant depending on the same param- 
eters as C3 that may vary from line to line. Proceeding as in the proof of Proposition 12.51 
the mass conservation and the energy equality fl3.8p imply the estimate fl3.12p and all 

bounds in fl3.10p - fl3.lip except those on dx\fh'e-, dxxhe, dx [he^^^^ and \/ed'^he. Using the 
previously obtained bounds (in particular the lower bound fl3.12p on he) and the interpo- 
lation inequalities 

\\dxxhe\\L°°{0,l) < C \\dl.he\\L'2{o,l) and \\dxhe\\L(i{Q^i) < C WdxxheW^l^Q^i^^ ll^2;^ellL2^(0,l) ' 
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one can estimate the last term on the right-hand side of the entropy equahty f l3.9p as 
follows: 



< C 6 sup 

[o,r] 



dxxUe ■ log he dx dt 



JQ 

T 



e'^\dxxU£\^ dx dt 







1/2 



dlK dxhedxJh ^ 
K hi 



sup 1 log heWl^io,!) 
[0,T] 

dxK\ 



hi 



L2(o,l) 



< 



Ce sup < II (9^/l£ 11^2(0,1) + II 11^2(0,1) ll^a^a^/l^l 1^00(0^1) + || Sa;/!^ H^e (q 1) f 
[0,T] J 



< Ce sup {||9^/le||L2(0,l) + C ||<9^/ie||L2(o,i) + C ||9^:./le||L2(o,i)} 
[0,T] 

+ C SUp||v/e5^/ie|U2(o,i) < CVi. 
[0,T] 



(3.13) 



Therefore, taking e sufficiently small, the bounds in f l3.10p -f l3.lip still to be proved follow 
from the entropy inequality (13.91) exactly as in the proof of Proposition 12.51 ■ 
We now show that solutions to the system fl3.4ap - fl3.4bp with boundary and initial 
conditions fl3.5p -f lXB]) converge to a solution to fl3.ip - fl3.2p as e — )■ 0. 

Theorem 3.3. For any positive a, Re,/3 and initial data {hQ,UQ) satisfying Ii2.9\) . there 
exists a global weak solution {h, u) to the system M.2a\) - ^1.2h\) with boundary conditions 



1.4\) - (T7B^) and initial conditions {ho,uo) in the sense of Definition \3.1i 



Proof: Take a sequence {£n}n>i ~^ and, for each n > 1, denote the corresponding 
solution to the approximate system fl3.4ap - p.4bp - fl3.5p - fl3.6p with e = by {he^.u^^). 

We investigate the compactness properties of the sequence {he„,U£^) and to this end 
use the estimates derived in Proposition 13. 2[ Let T > 0. First, owing to f l3.10p and 
f l332|) . i^/hZ) and {^/hZueJ are bounded in L°°((0,1) x (0,T)) and L~(0, T; ^^(0, 1)), 
respectively. Hence, by f]3.4bp . {dths„) is bounded in L°°(0, T; if~^(0, 1)). Since (/?.£„) is 
bounded in L~(0,T; H\0, 1)) and L'^{0,T; H\0, 1)) by flXTO]) and flXTT]) . it follows from 
the compactness of the embedding of H^{0, 1) in C([0, 1]) and Corollary 4 in Simon [lij 
that there is /i G C([0, 1] x [0,T]) such that, after possibly extracting a subsequence. 



he. 

dthe. 



^ h in L2(0,T;iy^'P(0,l))nC([0, 1] X [0,T]) for J9G[1,oo) 
^ dth in L°°(0,T;i7-i(0, 1)). 



(3.14) 



Combining f l3.12p and fl3.14p we additionally obtain that h is positive and 

Kl^h-^ in C([0,l]x[0,r]). 



(3.15) 
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We next turn to compactness properties of For that purpose, we write f l3.4ap as 

- ^ + BnKXhe. - eldtus„ (3.16) 

and claim that the right-hand side of (13.161) is bounded in L^(0, T; H^^{0, 1)). Indeed, by 
(13I0])-(13I2]), (Kyj, (v/^), and (v^^.m.J are in L-(0, T; ^^(0, 1)), L-((0,1) x 
(0, T)), and ^^^((0, 1) x (0,r)), which imply that {he^ulj and {dx{K„dxUeJ) are bounded 
in L2(0, T; H~^{0, 1)). Next, for any ip G C~((0, 1) x (0, T)) one obtains, using integration 
by parts and fl3^ - fl3:T2|) . that 



T rl 




Jo 

T 



i'he^dlhe^dxdt = d^xhe„[he„d^ip + ipd^hsj dxdt 

Jo Jo 

heJl^O,!) Ol^enl|L-(0,l)||5a-^||L2(0,l) + ||^||l-(0,1) ll^^^^n lU^CO,!)] 



< 



< c 



T 



m{o,i) 



1/2 



and 



nip he„dxll{he„) dx dt = dxiplli{hsjdxdt 
Jo Jo 

< ||ni(/lsJ||L-({0,l)x(0,T)) f / 11^11^1(0,1)^^ 



where Ili{h) is defined in (13. 3p . Similarly, 

i^K^dlhe^ dxdt 




'0 JO 
fT rl 



< e 



\\d%AL^{o,i) {U\\L'^(o,i)\\dlheAv\o,i) + \\heAL'^{o,i)\\dlip\\mo,i) 

,T;H'^(0,1))- 

Jo 

Finally, (m^J and {SndxxUe,,) are bounded in ^^((0,1) x (0,T)) and ^^(0, T; /f-^j^O, 1)), 
respectively, by (13. lip . Collecting the above information completes the proof of the bound- 
edness of the right-hand side of (I3.16p . whence 



{dtihe^Ue:)) is bounded in L\Q,T-H-\Q,1)). 



(3.17) 
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Now, owing to f l3.1ip and fl3.12p . we realize that is bounded in L^(0, T; if^(0, 1)). 

This fact along with fl3.10p - fl3.12p allows us to conclude that {h^^u^^) is bounded in 
L^(0, T; if^(0, 1)). Combining this with fl3.17p and Corollary 4 in Simon [lij ensures 
that {hs„Us„) is compact in L^((0, 1) x (0,T)). Then, thanks to f l3.15p . there exists 
u e L2(0, T; H\0, 1)) such that 

d^Ue^^d^u in L^{{0,1) x {0,T)) and Ue,, ^ u in L^{{0,1) x {0,T)). (3.18) 

The convergences fl3.14p . f l3.15p . and f l3.18p then allow to pass to the limit as n — )■ oo in 
the weak formulation of the approximating systems fl3.4ap - fl3.4bp - fl3.5p -f l3.6p in order to 
get that {h,u) satisfies (^-(^. ■ 

Remark 3.4 (Energy and entropy inequalities). It follows from the proof of Theorem 13.31 
together with the energy and entropy identities fl3.8p . f l3.9p . the property fl3.7p and the esti- 
mate fl3.13p that the weak solutions to fl3.ip -f l3^ obtained above satisfy the corresponding 
energy inequality 



sup 
te(o,T] Jo 
"1 



u 



Reh^ + U{h) + a 



dx + 



Jo 



< 



Re/^of + f/(M + ^^ 



dx 



u 



dxdt 
(3.19) 



and entropy inequality 



sup 

tG(0,T] Jo 



+ 



< 



-h{Reu + d^Sp{h)Y - -(f{h) + Re [a 
2 fj \ z 

i/io (Rewo + 5xV5(/io))^ - ^V5(^o) + Re (a^^-^^ 



+ U{h) 



dx 



dx. (3.20) 



Note that the right-hand side of fl3.19p - fl3.20p does not depend on T. Consequently, the 
statements of Lemma [231 and Proposition 12.51 hold also for the weak solutions to fl3.1l) - fl3.2p 
constructed in Theorem 13.31 



4 Limiting cases. 

4.1 Cases Re = and [5 = oo. 

By Remark 13. 4[ the statement of Proposition 12.51 is true for the weak solutions to (13.10 - 
(13. 2 p provided by Theorem 13.31 We may then investigate the behaviour of these solutions 
as either Re — ?■ or /3 — )■ oo. Let us first consider a sequence of positive real numbers 
(/3n), Pn — ^ oo, and denote the corresponding solutions to fl3.ip - fl3.2p with a fixed Re > 
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and /3 = /3„ by (/i/3„, 'U/?^)- Though the estimate on {up^/Pn) is useless in that case, one 
still recovers the estimate of in L^(0,T; -ffo(0, 1)) as a consequence of f l2.1ip . fl2.12p . 
and the Poincare inequality. Arguing as in the proof of Theorem I3.3[ we conclude that, 
after possibly extracting a subsequence, converges towards a weak solution to 

the strong-slip model f ll.7ap -f ll.7bp describing the dynamics of free suspended films. 

If we now consider a sequence of positive real numbers (Re„), Re„ — )■ 0, and denote 
the corresponding weak solutions to f l3.ip -f lX^ with Re = Re„ and a fixed (3 G (0, oo] 
by (^Re„5 "i^Ren), we may again proceed as in the proof of Theorem 13.31 to show that, after 
possibly extracting a subsequence, (/iRe„, MRe„) converges towards a weak solution to the 
strong-slip model f ll.2ap - fll.2bp without inertial terms. There is however an important 
difference as the bound on {dt{ /iRe„MRe„)) in -^^(0,T; if^^(0, 1)) is no longer available and 
we only obtain the weaker conclusion 

MRe„ in L2(0,T;ffi(0,l)). 

Still, owing to the strong convergence of (/iro„)) this allows to pass to the limit as n — )■ oo 
in all the terms which depends linearly on MRe„, that is, all the terms involving MRe„ except 
Re„/iRe„WRe„- But the latter converges to zero as Re„ — )■ since (/iRe„) and (uRe„) are 
bounded in L°°((0,1) x (0,r)) and L2(0, T; i/o^(0, 1)), respectively, by fliTT|) - fl2^ . 

4.2 Case cr = 0. 

As already pointed out in Remark |2.6[ for weak solutions (/i, u) to fl3.ip - fl3.2p . the estimate 
on dxh in (12. lip and both bounds in f l2.12p depend on a and explode as cr —)■ 0. In order to 
investigate the limiting behaviour of weak solutions to fl3.ip - fl3.2p as a — )■ 0, we refine the 
a priori bounds of Proposition 12.51 in the next proposition so as to avoid their dependence 
on 0", although paying for it with an exponential growth with respect to time. 

Proposition 4.1. For fixed positive o". Re, /3,T and initial data satisfying \2. 9\) . let {h,u) 
be a weak solution to liS. ip -z TOj) . There exists a positive constant C4 > 1 depending only 
on T, a, u, Re, E{uo,ho), and ^^(mo, /^o) defined in 112.1^) (but not on a & (0,1)^ such 
that the following terms are bounded by C4 in respective norms 

Vh, d^Vh, h-^/\Vhu G L~(0,T;L2(0,1)), (4.1) 

d.ih-'/^), V^d^^h, Vhd^u, u e L\{0,1) X (0,T)), (4.2) 

and 

C-^ < h{x, t) <C4, {x, t) e (0, 1) X (0, T). (4.3) 

Proof: First, the L°°(0, T; ^^(0, l))-estimates on Vh, Vhu and /i'^/^ ^nd the ^^((0, 1) x 
(0, T))-estimates on y/hdxU and u readily follow from the mass conservation and the en- 
ergy inequahty f l3.19p as before. Next, we actually estimate in a different way the term 
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Aull' {h)\dxh\'^ dxdt in the entropy inequality ( I3.20p . More precisely, using the estimate 
, Aa — 3h /, \ /, \ f 3 \ 

n(^) = ^5 > X{0,a){h)-^ + X{a,oo){h) j 

a 4 

where xa denotes the characteristic function of a set A C M and recalling that ^{h) = 
4i/ log h, one obtains 

uxkMki^ > - x,<.,»,(/.)»^ > - »^ . (4.4) 



Next, let us consider as in the proof of Proposition 12.51 the function Sp{u,h) defined in 
f l2.13p . Thanks to the mass conservation (11 .Op . the entropy inequahty f l3.20p also reads 

Sp{u{t),h{t)) + 



J J + Aau\d,j:,j:h\'^^ dx ds 

< Sp{uo,ho)-Au [ [ U'{h)\dxh\'^dxds. (4.5) 
Jo Jo 

Recalling that 

Sp{u,h)>\f h{Reu + dMh)fdx-^ (4.6) 
by ([2II5D, it follows from (jO]), (H^]), and Lemma [23] that 

-4z/ f Y\!{h)\dxh\^dx<-Aau [ dx + ^ [ h\dMh)\^ dx 

Jo Jo Jo 

< -Aau [ ^-^^dx+^ (si3{u,h) + ReE{u,h) + ^^ 

< -Aau + ^ ('^^^''' + ^^""O' + ^) ' 

where we used the energy inequality (I3.19P in the last estimate. Inserting this in (14. 5 p one 
ends up with 

/"* ( \d 
Si3{u{t),h{t)) + J J \Re— + Aau\dccxh\^ + Aau^j-^ \ dx ds 



4Re / 1 \ '± I 

< S[s{uo,ho)-\ { E{uo,ho) + — ] t-\ -/ S[s{u{s),h{s)) ds. 



17 



Applying Gronwall's inequality gives 



S^{u{t),h{t)) + I^Re^ + 4az/|(9^^/i|2^ dxdt 

< i Si3{uo,hQ) + Re 



E{uo,ho) + 



exp ( ) . (4.7) 



Combining (EZ]), (USD, and (iZl) implies the estimates on d^y/h in L°°(0,T; ^^(0, 1)) and 
dx{h-^/^), ^fad^Ji in L2((0, 1) x (0,T)) and completes the proof of ^^-^I^. To check 
(14. 3p . we first notice that the bound on ^JTi in L°°(0,T; -ff^(0, 1)) from (14. ip and the em- 
bedding of if^(0, 1) in L°°(0, 1) guarantee the upper bound in (14. 3p . Next, integrating the 
following identity 

1 1 r dMz) , , , 

with respect to y and using the Cauchy-Schwarz inequality give 

1 ^ dy ^ ( ndM^)\\ y f f dz 



K^) Jo Ky) \Jo Kz) J \Jo h^iz 

< I|/^"'/1K(o,i)+4||5xV^||l^(o,i)||/^"'/1l2{o,i). 

The lower bound in (14. 3 p then follows from the above inequality by (14. ip . ■ 
Using the a priori bounds from Proposition 14. H one can show existence of weak solutions 
to (I3l])-([321) in the case a = 0. 

Theorem 4.2. Let a = 0. For any positive Re,/3 and initial data {Jiq^uq) satisfying ^2. y\) 
there exists a global weak solution {h, u) to M.2a\) - ^1.2h\) satisfying \3. l\) -[37', 



Proof: Let us take a sequence {cr„} — t- and, for each n > 1, denote the corresponding 
weak solution to (I3.ip - (l3.2p with cr = (T„ by {ha„,Ua„)- Owing to Proposition 14. II we may 
proceed as in the proof of Theorem 13. 3 1 to show the existence of functions h and u satisfying 

h^^ h and h'^ h'^ in C{[0,1] x [0,r]), 
dtK„ ^ dth in L^{0,T;H-\0,1)) 

and 

d^u^^^d^u in L2((0,1) X (0,T)) and u^^ ^ u in L^{{0,1) x {0,T)). 
In addition, 

cTndx.K^^O in L2((0,l) x (0,T)) 
by (14. 2p . from which we deduce that 

(^ndxha„dxxha„ and (JuK^d^xKr, ^ ^ in L^((0, 1) x (0, T)) 

with the help of (I4.ip . This allows us to pass to the limit as ri — oo in all integrals in 
(ISII])-(E2]) and conclude that {h,u) satisfies ([XI])-([S2]) with a = 0. ■ 
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5 Convergence to solutions of intermediate-slip 
equation. 

In the introduction we mentioned that performing the change of variables (11 .Op and passing 
formally to the limit as /3 — )■ in f ll.2ap - fll.2bl) gives fll.lOa|) - fll.lObp and consequently that 
solutions of the strong-slip equation are expected to converge to that of f ll.Sp . The next 
theorem makes this formal procedure rigorous. More precisely, we show that scaled weak 
solutions to ffT:2aP - ffL2b|) converge to that of fll.Sp when a > 0. In this section, we denote 
the unsealed variables in f ll.2ap - fll.2bp with bars. 

Theorem 5.1. For fixed positive Re, a, let be a sequence of positive real numbers, 
Pn — )■ 0, and, for n > 1, denote a global weak solution to h3.1\) - [J7^) with initial data 
{ho,uo) satisfying /i2.9\) and (3 = (3n by (hn,Un)- Using U.9\) . we define 



K{x,t) ■= hnix 



/3 



and Un{x,t) := —Un\x 



(3 



G (0,1) X (0,oo). (5.1) 



Then there exist a positive function h and a subsequence of {hn,Un) (not relabeled) such 
that, for any T > 0, 

hn^h in L2(0,T;/7i(0,l))nC([0,l] X [0,T]), 
Un^u:=hd,{ad^^h-U{h)) m ^^((0, 1) x (0, T)), 

and h is a global weak solution to U.8\} considered with initial condition ho satisfying ^2. 9\) 
and boundary conditions ([i.y^[j-pT3]). 



Proof: Owing to f l3.ip -f lX^ . for each n > 1, the functions {hn,Un) satisfy the following 
scaled weak formulation 



oo /•! 



Jo 



hndtijj dxdt + / hoipi-, 0) dx 







oo fl 



Jo 



hnUnd^ip dxdt, 



(5.2a) 



oo fl 



JO 

oo fl 



hnUndtcj) dx (it + Re / hoUo(f){., 0) dx + Re 



oo fl 



JO 



hnu1^dx4> dx dt 



h„drU„dr(i> dxdt — (J 



JO 
oo fl 



JO 



JO 

{ahndxxhn - ^i{hn)) dr,(j)dx dt - 



dxhndxxhn4' dx dt 

oo ^1 

u^d dxdt = 



JO 



(5.2b) 



for all i) E C^i[0, 1] x [0, oo)) and G C^dO, 1) x [0, oo)). In addition, for any T > 0, 
{hn,Un) satisfies the scaled energy inequality 



1 r 



sup 

tG[0,T] Jo 
< 



U 



/3^Re/i„^ + f/(/i„) + a 



\dxhn\ 



dx + 



T fl 



JO 



[AuPn hn\dxUn\^ + ul) dx dt 



/3>/^o| + t/(M + ^^ 



dx 
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and the scaled entropy inequality 



sup / 

t6[o,r] Jo 

r-T rl 

+ 



(3: 



\dxhn\ 



+ U{K) - ^{h) 



dx 



^0 



{^nReul + Aav\dxxhn\^ + 4z/n'(/i„)|(9^/i„|2) dxdt 



</3n 



-/io(/3nRe uq + dx(p{ho)) + Re a 



}dccho\ 



dx — ^(ho) dx. 



Using the fact that a > and proceeding as in the proof of Proposition I2.5[ it follows from 
the mass conservation and the scaled energy inequality that 

(5.3) 
(5.4) 



{hn^^^)Adxhn) are bouuded in L°°(0, T; ^^(0, 1)), 
(un), {\/Ki\/KidxUn) are bounded in L^((0, 1) x (0,r)), 

and there exists C5 > 1 such that 

C^^ < hnix, t)<C5, (x, t) e (0, 1) X (0, T). 



We next use the scaled entropy inequality as at the end of the proof of Proposition 12.51 to 
estabhsh that 



{dx{h-^/^)), (dxxhn), are bounded in L\{0,1) x (0,T)). 



(5.5) 



Proceeding as in the proof of Theorems 13.31 and 14.21 one obtains that there exist functions 
h and u such that 

hn^h in L\0,T;H\0,l))nC{[0,l]x [0,T]), 
h-'^h~' in C([0,1] X [o,r]), 
u„^u in L\{0,1) X (0,T)), 
PnhndxUn^O ill ^^((0,1) X {0,T)). 

It is then straightforward to pass to the limit as /3„ — ?■ in fl5.2ap -( l5.2bp to obtain 



OO fl 



JO 

a 



hdt'4'dxdt+ / hQip{.,0)dx 



OO pi 



^0 



OO /"l 



^0 



dxxhdx{h(f)) dxdt + 



OO /"l 



^0 



hudxip dx dt, 

f OO pi 



Ili{h)dx4> dxdt 



~'0 



(5.6a) 

u(j)dxdt. (5.6b) 



Using an approximation argument, we may actually take = hdxijj in fl5.6bj) with ip G 
C~((0, 1) X [0, oo)) and obtain 



OO pi 



^0 



hdtip dx dt + / hQ^(.,0)dx 



a 



+ 



OO pi 



Jo 

OO pi 



Jo 



dxxhdx{h'^dx4') dx dt 
h'^dx(Il{h))ijdxdt, 
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which coincides with the weak formulation for f ll.Sp introduced in Bernis and Friedman 



[15] and Bertozzi and Pugh 19 



Remark 5.2. Owing to the positivity of the weak solution h to fll.8p constructed in 
Theorem I5.H equation f ll.Sp is then uniformly parabolic and it is likely that classical 
parabolic regularity results ensure that h is actually a classical solution to (11. 8p . As it is 



shown in Bertozzi et al. 17|] that there is a unique positive classical solution to (11. Sp . this 



implies the convergence of the whole sequence towards it. 



6 Discussion. 



We would like to note some open questions that remain to be solved for the strong-slip 
system (ll.2ap - fll.2bp as well as the corresponding shallow- water equations (ll.llap - fll.llbp . 
First, it would be of interest to figure out whether the weak solutions constructed in the 
previous sections are more regular. Partial answers have already been given for particular 
cases of (ll.llap - lll.llbp with a = 0, see Solonnikov ll| and Mellet and Vasseur |l3| where 
the existence of strong solutions is established. Furthermore, in Vaynblat et al. [20| the 
two-dimensional analogue of (II. Tap -i ll. 7bp describing dewetting of three-dimensional free 
suspended films is derived and reads 



dt{hu) + div {hu ®u) = div {h{Wu + (Vm)^)) 
dth = — div (hu) , 



2V(MivM) + cTWA/i- VP(/i) (6.1a) 

(6.1b) 



which corresponds to the choice fi{h) = h and X{h) = 2h in (ll.llap - fll.llbp . As mentioned 
in the introduction the BD entropy equality has been so far only established for (ll.llap - 
(ll.llbp if the relation (I1.13P for the viscosities is satisfied. Unfortunately this relation is 
not fulfilled for (I6.1ap - fl6.1bp . Analogously a two-dimensional version of (ll.2ap -f ll.2bp can 
be constructed but will present the same difficulty, so that the analysis presented herein is 
not likely to extend in a straightforward way to the two-dimensional model. 
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